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The early idea that a non-perturbative gluon condensate affects the spectrum of heavy quarks is revisited in 
the light of modern simulation techniques. We evaluate the low lying spectrum of bound states of two heavy 
quarks for large hypothetical quark mass, ttiq > mj, using non-relativistic QCD and compare with other models 
to test the consistency. 



The concept of a non-perturbative gluon conden- 
sate was initially introduced by Shifman, Vain- 
stein, and Zakharov Q. A non- vanishing gluon 
condensate can give rise to a shift of energy lev- 
els of heavy quarkonia by second order Stark 
splitting. Leutwyler and Voloshin [|j calculated 
the level shifts in leading order perturbation the- 
ory. Recent improvements like the stochastic 
vacuum model SVM [|, non-perturbative field 
strength correlators |3J or a combination of the 
SVM with lattice calculations try to solve 

some of the problems of the original LV-model. 
Using NRQCD has the advantage that no apri- 
ori knowledge about the heavy quark potential 
or the vacuum background field is needed. The 
physical problem is reduced to a simulation of 
non-relativistic heavy quarks in a correlated, rel- 
ativistic gluonic background. We use an 0(a 4 ) 
improved gauge action with a plaquette and a 
1x2 Wilson loop term for the relativistic simula- 
tions and an 0(v 4 , v 6 ) corrected non-relativistic 
action 

S NR = ^ (D t + H +SH(c t )) 

with tree level couplings, Cj = 1, for the heavy 
quarks. The quark's Green's functions are cal- 
culated from the standard (asymmetric in 5H) 
evolution equation. Meson propagators with spe- 
cific quantum numbers are obtained from the 
quark propagators using a smearing procedure 
with trial wave- functions from a solution of the 
radial Schrodinger equation using a method from 
ref. Jl3| for each value of the bare quark mass. 
Preliminary results based on Coulomb wave- 
functions and a small subset of the configurations 
were already reported at a previous conference 



For the simulations we use a lattice of size 16 4 
and a gauge coupling of /3 = 9.17 to obtain a 
very fine spacing necessary to resolve small bound 
states. After generating 360 decorrelated configu- 
rations on a Quadrics Q4 we fix to Coulomb gauge 
with a stopping criterion, 9 < 10~ 6 . The simu- 
lations require 60 days of Quadrics Q4 time, 70 
days on a workstation for gauge-fixing, and more 
than 150 days of CM2/CM5 time. We apply stan- 
dard correlated multi-state multi-exponential fits 
with vector fits for a set of S states and ratio 
fits for hyperfine and S-P splittings. The error 
is determined by a jackknife procedure including 
the uncertainty from the fit range, [t m i n , N T ]. 
A consistent scale is obtained from V q q and the 
force indicating that scaling violations are indeed 
small at this high value of the gauge coupling (see 
ref. for details). 

As a first step our results can be compared to 
leading order perturbation theory which allows to 
separate the Coulomb from the condensate con- 
tribution. Following refs. HQ the level splittings 
are given to lowest order by the expressions 
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There are two length scales in the system: the 
gluon correlation length Xq and the quark corre- 
lation length Xq which is related to the rotational 
period of the two heavy quarks inside the quarko- 
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Figure 1. The hyperfine splitting between the 
ground state of the 3 Si and the 1 S'o meson. The 
line gives the perturbative Coulomb contribution 
with aJrriQ) as shown in fig. 3. Also shown are 
the cc flllll and bb values |[12fl. 



nium, Aq ~ 27r^/ (r 2 )/(u 2 ). The shortest correla- 
tion length governs the behaviour of the energy 
levels. There are two limiting cases: a static con- 
densate (Ag >> Aq) when f n i = 1 and a rapidly 
varying condensate (Aq << Aq) when /„; < 1. 

We can solve the set of equations {AEhf s , 
AEsp} for the two unknowns a s and (a s G 2 ) in- 
serting the values for the splittings from NRQCD 
shown in figs. 1 and 2. In this way we obtain the 
strong coupling and the condensate as a function 
of the bare quark mass shown in figs. 3 and 4. If 
we reverse this procedure for a while and accept 
that fig. 3 describes a reasonable behaviour of the 
running coupling a s (m(j) then the perturbative 
Coulomb contribution (~ 1.05 mQaf) shown in 
fig. 1 as a line together with our data completely 
dominates the hyperfine splitting Ehf s . On the 
other hand for the S-P splitting the condensate 
contribution is larger than the Coulombic part 
which has been criticised before and sheds doubt 
on the treatment of the condensate as a pertur- 
bation. The magnitude of the contributions to 
Ehfs and E$p does not change significantly as 
the gluon condensation length is varied. The 
parameter \q merely rescales the value of the 
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Figure 2. The S-P splitting from NRQCD 
(squares), perturbation theory and for the 
Buchmiiller and Tye model Q . 

condensate by a factor /^(Ag). If we neverthe- 
less take the perturbative results seriously we see 
from fig. 4 that it is possible to get a value of 
the condensate compatible with sum rule results. 
However, the value of (a s G 2 ) strongly varies with 
quark mass. The static approximation (Ag — » oo) 
should be valid for large quark mass while the 
approximatiuon of a rapidly varying condensate 
(Aq << Aq) is expected to hold for intermediate 
masses, tuq < 40 GeV, because Aq is decreasing 
with rriQ. 

Summary 

At large tuq the hyperfine splitting of the 
ground state becomes compatible with Coulom- 
bic behaviour. However, the S-P splitting from 
NRQCD is neither in accord with the B+T model 
nor with leading order perturbative results for 
LV-type models. On the level of lowest or- 
der perturbation theory no consistent picture of 
the effect of the gluon condensate can be found 
allthough there is evidence that the general trend 
predicted by the analytic results is qualitatively 
correct. The NRQCD hyperfine splitting of the 
ground state is decreasing and the S-P splitting 
is increasing with tuq for large quark mass con- 
sistent with the perturbative dependence on dif- 
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Figure 3. The running coupling a s from a com- 
parison of the S-P and the hyperfine splitting of 
the ground state from NRQCD with perturbative 
predictions. 

ferent powers of a s . It will be interesting to see 
15 1 if a comparison with more rigorous models 
like the combined SVM and lattice approach will 
consolidate the situation. An alternative might 
be a closed lattice approach. Then one is faced 
with the difficult problem to switch the conden- 
sate on and off during the simulations. 

Acknowledgement 

J.F. was supported by a DFG fellowship during 
the initial phase of the project. We also wish to 
thank the IAI at the University of Wuppertal and 
the computer centre of the University of Erlangen 
for their support. 

REFERENCES 

1. V.I. Zakharov, (1978); M.A. Shifman, 
A.I. Vainstein and V.I. Zakharov, Nucl. Phys. 
B147 (1979) 385 and 448. 

2. M.B. Voloshin, Nucl. Phys. B154 (1979) 365; 
H. Leutwyler, Phys. Lett. B98 (1981) 447. 

3. A. Di Giacomo and H. Panagopoulos, Phys. 
Lett. B285 (1992) 133. 

4. W. Buchmiiller and S.-H.H. Tye, Phys. Rev. 



0.2 



0.15 



0.1 



0.05 



; G 2 > [GeV 4 ] ' 1A G =0'^' 
1/?L G =350 MeV — x— 

svz 

Narrison 




20 25 30 35 40 45 50 
m Q [GeV] 

Figure 4. A comparison of the gluon condensate 
for the static case, Xq — * oo and for A^ 1 = 350 
McV with different values from the literature 
(SVZ: |l|, Narrison: @). 
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